Structure formation with a self-tuning scalar field 
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A scalar field with an exponential potential has the particular property that it is attracted into a 
solution in which its energy scales as the dominant component (radiation or matter) of the Universe, 
contributing a fixed fraction of the total energy density. We study the growth of perturbations in 
a CDM dominated Q = 1 universe with this extra field, with an initial flat spectrum of adiabatic 
fluctuations. The observational constraints from structure formation are satisfied as well, or better, 
than in other models, with a contribution to the energy density from the scalar field fi^ ~ 0.1 which 
is small enough to be consistent with entry into the attractor prior to nucleosynthesis. 



The simplest viable cosmology which follows from in- 
flation, a flat universe with pressureless matter and 5% 
baryonic dark matter, has been unable to fit both the cos- 
mic background radiation (CBR) fluctuations and mea- 
surements of mass fluctuations on scales of a few Mega- 
parsecs. The paradigm of inflation is sufficiently com- 
pelling that there have been various attempts at modi- 
fying this 'standard cold dark matter' (sCDM) scenario 
[jij. The possibility that some part of the energy den- 
sity of the Universe is in a form other than particle-like 
matter has been envisaged, in particular in the form of 
a constant energy (ACDM) j| or time-dependent coher- 
ent energy density in a scalar field |^^]. In this letter 
we discuss the cosmology of a model with a scalar field 
which has a simple exponential potential. It is distinctly 
different from other scalar field cosmologies, in that its 
energy density plays a role from very early times, rather 
than just at recent epochs, and resembles much more the 
'mixed dark matter' (MDM) model Q in which there is 
a component of matter which is collisionless during a pe- 
riod of the growth of structure. The required potential 
has the merit that it arises quite generically in particle 
theories involving compactifications such as supergravity 
or superstring theories, and has (mainly for this reason) 
been quite extensively discussed in the context of infla- 
tionary models. 

Let us first explain the properties of an exponential po- 
tential which make it a particular and interesting case. 
The equations of motion in an expanding FRW universe 
for the homogeneous mode of a scalar field <f> with po- 
tential V{4>) coupled to ordinary matter only through 
gravity are 
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where p n is the energy density in radiation (n = 4) or 
non-relativistic matter (n = 3), 7i = - is the conformal 
expansion rate of the universe with scale factor a, dots 
are derivatives w.r.t. conformal time r, ' = 4r and Mp = 



2.4 x 10 18 GeV is the reduced Planck mass. Multiplying 
(0) by 4> and integrating, one obtains 



;{a) = p{a )e J »° 



(4) 



where P<h = 



2 + V(4>) is the total scalar energy, and 
£ = V(<p)/ ' p$. In general therefore the energy density of 
a scalar field has the range of possible scaling behaviours 
p oc 1/a™ 1 with < m < 6, and the scaling is completely 
determined by the ratio of its potential to its kinetic en- 
ergy. 

The special cosmological solutions in which we are 
interested here are attractor solutions of (|l|) - (fl) for 
the case of an exponential potential V(<f>) = V e~ x ^/ Mp , 
which were given in || and f|. In these solutions the 
scalar field evolves so that its total energy density p$ 
scales in the same way as the dominant component (i.e. 
p$ cx 1/a") and contributes a fixed fraction of the total 
energy density given by 
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for A > 1/ yfn. Note that it is A alone which determines 
the solution. The existence of the attractor can be un- 
derstood to follow from the fulfillment of two conditions: 
(i) scales faster than l/a n if p^ » p n and, (ii) scales 
slower than 1/a™ if p^ << p n . These two behaviours 
tend to drive the two components to the attractor which 
lies between them. That the first condition is satisfied 
can be seen from solving (Q) - (|J) with p n = for the 
exponential potential. There is then a different set of 
attractors B in which 
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where A < -\/6-For A > \/& there is not a single attrac- 
tor, but all solutions have £ —> asymptotically (and, 
therefore, p oc I /a 6 ). The condition A > l/y/n for the 
attractor is indeed therefore just that anticipated. 
The second condition can be understood qualitatively as 
follows. Taking, for simplicity, the case n = 4, (|l]) - (|^) 
with V{4>) = give <fi(t) = o (^f ) 2 and therefore 
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These will also hold as approximate solutions in the case 
that the potential energy is sub-dominant. The first so- 
lution shows how, for a sufficiently steep exponential, the 
potential energy can remain small relative to the kinetic 
energy (~ 1/t 3 ) so that the rapid scaling (associated 
with £ << f) can be maintained. The field in both the 
attractor solutions (@) and (||) has this same logarithmic 
time dependence. On the other hand, the second limit 
shows how the larger damping due to radiation domina- 
tion slows down the evolution of the field giving an al- 
most constant potential energy which will thus ultimately 
catch up with the kinetic energy, increasing £ and causing 
the scalar energy to scale slower. 

A potential which is less steep than the exponential 
will not satisfy the first condition 0, and a steeper po- 
tential (e.g. ~ I m p) will always decay asymptot- 
ically relative to the other components. The existence 
of this particular attractor cosmological solution is thus 
quite specific to the exponential potential. Further this 
is in fact a potential which can arise quite generically in 
particle physics theories involving compactified dimen- 
sions (with internal dimensions characterized by Mp). 
For this reason it has been considered quite extensively 
in the context of inflation |[8|— p~0| , since for A < y/2 the so- 
lutions (||) describe 'power-law' inflation (with a oc t 2 l x 
in terms of physical time t = J adr) . Examples of spe- 
cific supergravity theories in which such potentials are 
obtained are given in , and various higher dimensional 
theories of gravity in which they arise discussed in detail 
in 0] and 

If such a field does exist, it will enter the attractor 
and contribute a fraction of the energy density (fixed by 
A) at some time determined by its initial energy density. 
Nucleosynthesis provides the earliest constraint on how 
large such a contribution can be. The expansion rate 
of the Universe at nucleosynthesis is increased over its 
standard model value by the same amount as A N e f f rel- 
ativistic degrees of freedom, with fi</, = 1 10 J^^&N^f/A » 
where D,^ is the fraction contributed in the matter era. 
There is some disagreement on the precise nucleosynthe- 
sis constraint on AN e ff, but a bound of AN e ff = 0.9 is 
given by various authors |l2j or even a more conservative 
one of AN e ff = 1.5 by others |L3]], which corresponds to 
fi < 0.1 - 0.15. 

Prima facie this constraint would seem to require en- 
try into the attractor after nucleosynthesis if the scalar 
field is to play any significant role cosmologically ||] . The 
requirement of entry after nucleosynthesis would appar- 
ently mandate the unattractive fine-tuning (typical of 
scalar field models) of the initial energy density in the 
potential to some small value. It was in fact the incor- 
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FIG. 1. Mass variance per unit Ink computed from Boltz- 
man code for different models compared with that inferred 
from a compilation of galaxy surveys |l9| 



rectness of this second assumption which motivated the 
present study: If, prior to nucleosynthesis, the energy 
density in the exponential field with A > V6 dominates 
over that in the radiation, there will typically be a long 
transient period after ~ p 7 during which the scalar en- 
ergy is very sub-dominant (much less than its value in the 
attractor (||)). This is simply because the ratio £ — > in 
the kinetic energy dominated pure scalar cosmology, but 
is of order one in the attractor with radiation. During 
the time that £ is increasing (potentially many expan- 
sion times as it cannot grow faster than a 6 ) the scalar 
field energy continues to red-shift away as 1/a 6 . Such a 
dominance by kinetic energy can occur in certain post- 
inflationary cosmologies which have considerable interest 
in their own right |14[] , [jfJJ. It has transpired from the 
present work however that the first reason for disregard- 
ing this model is also incorrect, and that entry to the 
attractor prior to nucleosynthesis is in fact consistent - 
quite simply because the small contribution has a com- 
pensating long time to act. 

We have carried out a detailed calculation of the evolu- 
tion of perturbations in this cosmology (which we refer to 
as <^CDM). We assume that the attractor is established 
at the beginning of our numerical simulation, deep in the 
radiation era, and take an initial standard inflationary 
scale-invariant spectrum of adiabatic perturbations. The 
relevant equations are the linearized coupled Einstein- 
Boltzmann equations given in |l6f , supplemented by the 
scalar field and its perturbations 4>totai — 4>{ T ) + <p( T > x )> 
with evolution equation 
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and additional components to the perturbed energy- 
momentum tensor: 
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where 7 is the trace of the metric perturbation. We 
vary fi^ and h (where Ho — /il00km/s/Mpc is the Hub- 
ble constant today), keeping the remaining cosmologi- 
cal parameters fixed at the values of sCDM, and find 
the best fit model to both CMB and large scale struc- 
ture. To do this we use the COBE measurement of CMB 
anisotropics on large scales [ p"7| to normalize our theory 
p8[ , estimate the theoretical mass variance per unit Ink, 
Ar(k) and compare with that rendered from a collec- 
tion of galaxy surveys |19|. In Figure [l] we show A 2 (k) 
for two best fit 0CDM models, for sCDM, for a ACDM 
universe with fl\ = 0.6, and for an MDM model with 
f2„ = 0.2 in the form of two massive neutrinos species. 
It is clear that for these values 0CDM fares as well or 
better than the other models. Another useful quantity 
to work with is the mass fluctuations on 8h _1 Mpc scales, 
a \ = _£» fA 2 {k)(^^) 2 U =8 . This can be related to 
masses and abundances of rich clusters and supplies us 
with a very tight constraint on possible cosmologies; in- 
deed current estimates give ag = 0.6 ± 0.1 p0| . A good 
fit to as is 
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where af DM is the COBE normalized sCDM erg. Again 
we see that there is range of values of £7^ and Hq which 
satisfy the above constraint and are consistent with the 
limits imposed by BBN. In Figure || we compare the Ces 
of our models with a compilation of data points |^T| . 
Again they are consistent with the current data. 

The evolution of perturbations in the presence of the 
scalar field is simple to understand. On superhorizon 
scales there is the usual growing mode with 6 C) tp oc t 2 
(where S c is the density contrast in the CDM). This is 
to be expected; the superhorizon evolution is insensitive 
to the "chemistry" of the matter and totally dominated 
by gravity. On sub-horizon scales in the radiation era, 
the Meszaros effect comes into play giving S c oc lnr. The 
specific effect of the scalar field appears on subhorizon 
scales in the matter era. The perturbation in the scalar 
field itself has the approximate solution <p oc -^j^ J 3 (fcr) 
(where J v is a Bessel function) which, when fed back into 
the equation for S c gives an altered solution for the usual 
growing mode 5 C oc r 2 ~ e where 
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This solution shows explicitly how even a small contribu- 
tion from the scalar field can give a significant effect, as 
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FIG. 2. Comparison of different model predictions to cur- 
rent experimental data. All models were COBE normalized 
and are labeled as in Figure 1. 



it acts all the way through the matter era. The expected 
suppression of \8 C \ 2 for modes larger than k eq is of order 
(1 + z eq )~ e , where k eq is the wavenumber of the hori- 
zon size at radiation-matter equality. This last effect is 
reminiscent of the evolution of perturbations in a mixed 
dark matter (MDM) universe where one has component 
of matter, p u which is collisionless for a period of time 
during the matter era [p2| . 

It is useful to pursue a comparison between 0CDM and 
MDM to identify the key differences. Firstly the scaling 
behaviour of the additional background energy density 
differs: While for </>CDM the energy density in <j> follows 
the dominant form of energy quite closely, for MDM p u 
changes from scaling as 1/a 4 to scaling as 1/a 3 when 
iksTu ~ rriy where T v (m u ) is the massive neutrino tem- 
perature (mass) and ks is the Boltzmann constant. For a 
period between matter-radiation equality and this transi- 
tion Q v is smaller than its asymptotic value, and there is 
less suppression of growth in the CDM than in the case 
of the scalar field. A further difference is that the pe- 
riod of time during which perturbations are suppressed is 
shorter in MDM compared to </>CDM. In both cases there 
is a wavenumber k su which separates growing modes from 
damped modes. For 0CDM this scale is roughly the hori- 
zon i.e. k su oc i, while for MDM it is the free streaming 
scale i.e. k su = 8a 1 1 2 (m v / 10eV)hMpc~ 1 oc r. Clearly in 
the latter case any given mode of 5 C will eventually start 
to grow. In particular modes around k eq will already have 
started to undergo collapse. A final important difference 
concerns the evolution of perturbations in the radiation 
era. For MDM, the perturbation in the massive neutrinos 
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behaves much like radiation until it is well inside the hori- 
zon, and this transition is set by the Jeans scale i.e. when 
hr ~ l/c s = For the scalar field, on the other hand, 
the transition occurs for larger wavelengths, fcr < 1. This 
means that perturbations in the CDM will stop growing 
earlier in 0CDM than in MDM. The accumulated effect of 
these differences explains what we have observed - that, 
with half the energy density of MDM with two massive 
neutrinos, </>CDM brings about approximately the same 
suppression of power on small scales. The fact that this 
suppression lasts until today leads to the formation of 
structure at higher redshifts. There is now strong obser- 
vational evidence that any model of structure formation 
must have this feature. 

Let us now turn to the effect that the scalar field has on 
the CMB. We shall rely on the simplified picture of p3| to 
understand the angular power spectrum, Ci, defined as 
C(9) = (^(n)^(n')) = (4^-^(2^+ l)C t P e (caa 6), 
with n ■ n' = cos 9. For t > 100 the main features of the 
CiS are given by the power spectrum of radiation pertur- 
bations at last scattering, <$ 7 . Ignoring projection effects, 
one has that the structure of the peaks and troughs are 
given by 
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where r s is the sound horizon in the baryon-photon fluid, 



dr 
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k is roughly related to the angular frequency I. The fact 
that the properties of the Cgs are dominated by this quan- 
tity at a ~ 10 -3 means that the effect of <fi on the CMB 
will be much smaller than its net effect on <5 C . Adding 
the scalar field component brings about two effects which 
we can understand qualitatively. Firstly the oscillations 
are shifted to higher is. Because of the additional en- 
ergy density in the scalar field, the expansion rate will 
be larger and the conformal horizon will be smaller for 
the same red-shift in 0CDM compared sCDM. This feeds 
through to give a different r s for the same value of a, 
shifting the peaks as observed. The other main feature 
is an increase in power in the peaks. This can be un- 
derstood easily using the picture outlined in |e3]. The 
oscillations in <L are driven by the evolution in the grav- 
itational potentials, and here as in the MDM case 63 ] 
the change in the growth of metric perturbations boosts 
the amplitude of the peaks by a few percent. 

We conclude that the cosmological model we have stud- 
ied provides an interesting and distinct alternative to 
other models which have been proposed. It has the at- 
tractive feature that A (= the single extra pa- 
rameter compared to standard CDM, has a value which 
is of the order naturally expected in the many particle 
physics theories in which the field arises. With the launch 
of high resolution space based experiments, such as the 
Planck explorer and the MAP satellite, it should be possi- 
ble to distinguish the effect on the CMB of such an expo- 
nential scalar field if it exists, or to rule out its existence 



and place tighter constraints on the physical theories in 
which these fields arise j25| . 
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